Low-momentum interaction in few-nucleon systems 
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The effective low- momentum interaction Viowfc is applied to three- and four-nucleon systems. 
We investigate the 3 H, 3 He and 4 He binding energies for a wide range of the momentum cutoffs. 
By construction, all low-energy two-body observables are cutoff-independent, and therefore, any 
cutoff dependence is due to missing three-body or higher-body forces. We argue that for reasonable 
cutoffs V| ow fc is similar to high-order interactions derived from chiral effective field theory. This 
motivates augmenting Viowfc by corresponding three-nucleon forces. The set of low-momentum two- 
and three-nucleon forces can be used in calculations of nuclear structure and reactions. 
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Microscopic nuclear many-body calculations are com- 
plicated by the short-distance repulsion in nuclear forces, 
which leads to strong high-momentum components in nu- 
clear wave functions. Usually, one solves this problem 
by introducing an effective interaction, the Brueckner G 
matrix, which resums in-medium particle-particle scat- 
tering. The G matrix is a soft interaction, which is both 
energy- and nucleus-dependent and typically requires ap- 
proximations in practice. Moreover, the resummation of 
particle-particle contributions makes it extremely com- 
plicated to treat particle-particle and particle-hole corre- 
lations on an equal footing. 

An alternative strategy to construct a soft interaction 
by integrating out the high-momentum components in 
free space has been formulated in Using a renor- 
malization group (RG) approach, phenomenological two- 
body potential models can be evolved to an effective low- 
momentum interaction, called Vi owfe , which is energy- 
independent, hermitian and preserves the on-shell T ma- 
trix below a cutoff A in momentum space as well as the 
deuteron binding energy. For A < 2fm~ 1 , the matrix 
elements of V\ ow k are practically independent of the po- 
tential model it is derived from and thus unifies all nu- 
clear forces used in microscopic nuclear structure calcula- 
tions Q. By construction, Mowfe is much softer than the 
modern potential models, and thus can be used directly 
for microscopic nuclear calculations in different mass re- 
gions 3, 4] or for different densities 0,0 This is clearly 
important to theoretically extrapolate to the nuclear drip 
lines without ambiguities due to unknown interactions. 

Over the last few years, there has also been an immense 
progress in our understanding of nuclear interactions 
from chiral effective field theory (EFT). This approach 
qualitatively explains the hierarchy of two-nucleon (2N), 
three-nucleon (3N) and higher-body forces 0, which is 
observed using phenomenological models. On a quanti- 
tative level, it was shown that the resulting 2N and con- 
sistent higher-body interactions lead to a qui te g ood de- 
scription of 2N as well as 3N observables |lBllSlIIllI3 
In the pionfull EFT approach, the Lippmann-Schwinger 



equation is regularized by imposing a cutoff A w 2.5 — 
3.0 fm -1 . Thus, the chiral potentials are also low- 
momentum interactions, and with the universal property 
of V[ ow k , this suggests that Vj ow k effectively parameter- 
izes higher-order chiral 2N interactions. While EFT of- 
fers the only systematic approach to consistent 2N and 
higher-body forces, Vj ow fc can be evolved to arbitrary cut- 
offs with cutoff-independent 2N observables. 

Since Viowfe is constructed within the 2N system, one 
neglects many-body forces due to degrees of freedom 
missing in the effective theory (contributions from the 
A) as well as due the truncation to low momenta (con- 
tributions from high- momentum nucleons). In any effec- 
tive theory, these effects are inseparable. In this Letter, 
we use cutoff dependence as a tool to assess the effects 
of many-body forces. Motivated by the similarities be- 
tween Viowfe and chiral low-momentum interactions, we 
combine Vi ow ft with the leading chiral 3N force to absorb 
the cutoff dependence in A ^ 4 binding energies. Fi- 
nally, we examine the expectation values of the various 
force components to check that the hierarchy of nuclear 
two- and three-body forces is maintained. 

We first calculate 3N and 4N binding energies by solv- 
ing the Faddeev-Yakubovsky equations with only the 
two-body Viow fc ■ We include electromagnetic and isospin- 
breaking effects and vary the cutoff over a wide range. 
Our results are numerically stable for the studied cutoff 
values, which requires a careful treatment of the neces- 
sary interpolations in the vicinity of the sharp cutoff. We 
also checked the convergence with respect to the included 
partial waves. We estimate an accuracy of 2 keV for the 
3 H and 3 He and 50keV for the 4 He calculations. More 
details about the numerical method can be found in |l3| . 

In Fig. ^ we give results for binding energies of the 
3N system. We show results for the Vi ow fc derived from 
the CD-Bonn 2000 Q and Argonne v w [3 interac- 
tions. The cutoff dependence is due to missing three- 
body forces. For large cutoffs, we reproduce the known 
binding energies obtained with the bare interactions only. 
For intermediate cutoffs, we find a stronger binding with 
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Mowfe- This could be expected, because softer interac- 
tions generally lead to stronger binding. It is also con- 
sistent with the correlation between the triton binding 
energy and the deuteron D-state probability observed 
for phenomenological potentials 16] . For Viowfc the D- 
state probability decreases monotonically with a decreas- 
ing cutoff. Therefore, this correlation evidently breaks 
down for cutoffs below A w 1.6 fm - . The binding then 
decreases, as attractive parts of the bare interactions are 
integrated out. 

For cutoffs A < 2m w , truly model- independent results 
are obtained and the binding energy curves for the CD- 
Bonn 2000 and Argonne Vi$ Vj ow fc interactions collapse. 
In Fig.^ w e also show the cutoff dependence of the differ- 
ence in 3 He and triton binding energies, which is due to 
electromagnetic and isospin-breaking contributions. The 
difference varies by 60 keV and correlates with the bind- 
ing energy, since the latter is related to the charge ra- 
dius For special choices of the cutoff, both exper- 
imental binding energies can be reproduced simultane- 
ously without a 3N interaction. We emphasize that 3N 
forces will contribute to other observables. Nevertheless, 
it may be interesting to study many-body systems using 
these particular values for the momentum cutoff, since a 
simple zero-range 3N force vanishes in these cases. 

Our results indicate that 3N forces due to the trun- 
cation to low momenta are of the same order as ad- 
justed 3N forces due to missing excitations of nucleons, 
although these effects cannot be separated. The bare 
3N forces provide about 0.7—1 MeV of binding in con- 
ventional models, whereas the binding energies given by 
Viow k change by 1 MeV over the large cutoff range. In 
this sense the truncation to low momenta does not in- 
duce strong three-body forces. We note that this is in 
contrast to the interpretation given in [l^. There, the 
size of 3N forces was assessed by comparing the Vj ow fc 
binding energies to the results of the bare 2N potential 
model. This neglects the uncertainty in the binding en- 
ergy predictions of traditional 2N forces and misses that, 
in effective theory approaches, the effects of the trunca- 
tion to small cutoffs are inseparable from those of missing 
degrees of freedom like the A. Because these two contri- 
butions to higher-body forces cannot be disentangled at 
low energies, we will absorb both by augmenting Vi ow fc 
with a chiral 3N force below. 

For further insight, we have calculated the a-particle 
binding energy. To obtain an overview, calculations 
are performed for the smallest cutoff considered A = 
l.Ofm - , in the maximum of the triton binding energy 
at A = 1.6 fm~ , for two cutoffs which lead to 3N bind- 
ing energies close to the experimental one, A = 1.3fm _1 
and A = 1.9 fm -1 (Argonne i>i$) or A = 2.1 far -1 
(CD-Bonn 2000), and for a cutoff in the tail at A = 
3.0 far 1 . The focus of our studies is whether the cutoff 
dependence of Viowfc can be related to correlations ob- 
served when traditional two-body interactions are used. 
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FIG. 1: (Color online) Cutoff dependence of the 3N binding 
energies and the binding energy difference of 3 H and 3 He. 
Results are shown for the Argonne Dig and the CD-Bonn 2000 
potential. The horizontal solid lines represent results for the 
bare two-body interactions and the dotted lines denote the 
experimental binding energies. 



From fliil flfl Eoj ] , it is well-known that there is an al- 
most linear relation between 3N and 4N binding ener- 
gies, known as the Tjon-line. This correlation holds with 
very good accuracy for all modern interactions, but is 
slightly broken by the action of 3N forces. As can be 
seen in Fig. the various V\ Q v/k results do not differ 
significantly more from the phenonemological Tjon-line 
than calculations with adjusted 3N forces. We see that 
as a further indication that 3N and 4N contributions are 
not unexpectedly large due to the low-momentum trun- 
cation, at least for the triton and a-particle. Already at 
A = 3.0 fm -1 the Vj wfc prediction is almost exactly on 
the Tjon-line given by the phenomenological models. 

As also seen in Fig. even if a cutoff is chosen that 
leads to a good description of the 3N binding energies, 
the 4N binding energy deviates from experiment. Clearly, 
3N or higher-body forces must act for these values of the 
cutoff. In the following, we construct a low-momentum 
3N interaction by fitting the leading chiral 3N force to 
Viowfc- For simplicity we restrict ourselves to the Viowfc 
derived from the Argonne v±g potential. The chiral 3N 
force to leading order contains a long-range 27r exchange 
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FIG. 2: (Color online) Correlation of the 3 H and 4 He bind- 
ing energies. Results are shown for several modern poten- 
tial models alone (plusses) and with adjusted 3N forces (di- 
amonds) [20|. The Viowk results are for the Argonne vis 
(squares) and the CD-Bonn 2000 potential (crosses). The 
solid line is a linear fit to the 2N force model results only. 



A [fm- 1 ] 


CD 


CE 


1.0 


3.621 


5.724 


1.3 


11.889 


2.265 


1.6 


2.080 


0.230 


1.9 


-1.225 


-0.405 


2.5(a) 


-0.560 


-0.707 


2.5(6) 


-3.794 


-1.085 


3.0(*) 


-7.500 


-2.151 



TABLE I: Fit results for cn and ce for various cutoffs 
of the Viowfe derived from the Argonne ms potential (for 
(*) see text). The strength of the 2-n exchange part is de- 
termined by ci = -0.76 GeV"\ c 3 = -4.78 GeV -1 and 
c 4 = 3.96 GeV -1 HJ. 
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FIG. 3: (Color online) Relation between Co and ce obtained 
by requiring that V\ wk augmented by the 3N force predicts 
the 3 H binding energy correctly. 
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FIG. 4: (Color online) Dependence of the eigenvalue r\ of the 
Yakubovksy equation on cd for various cutoffs. A deviation 
of 7} — 1 = 0.01 corresponds to a deviation of approximately 
600 keV from the experimental value. 



part, an intermediate range one tt exchange (D-term) and 
a zero-range contact interaction (£7- term), see 0, H^ . 
For the operator form and the definition of the strength 
constants, we refer the reader to Eqs. (2) and (10) in |12]| . 
The interaction is regularized by exponential cutoff func- 
tions of the form exp(— (j?/A) 8 ) with the cutoff taken 
from Viowfe- The very high exponent guarantees a very 
sharp drop to zero at p = A. The 2n exchange part 
is determined by strength constants Cj, which we take 
from 21], where they were obtained by a fit to NN 
data. The dimensionless strength constants cjj and Ce 
were obtained from a fit to the 3 H and 4 He binding ener- 
gies. First, a relation between cd and ce was established 
by requiring that the 3 He binding energy of —8.482 MeV 
is described accurately. The resulting dependence for 
various cutoffs is shown in Fig. |21 For small cutoffs we 
obtain a linear relationship, which suggests that the D- 
and E'-terms are perturbative in this region. We have 
checked explicitly and also for the c-terms that these are 



perturbative for A < 2fm~ . This could be useful for 
applications, where it is practically impossible to include 
the 3N force into the dynamical equations, but a pertur- 
bative treatment is feasible. 

In Fig. 21 we show the eigenvalue 77 of the Yakubovsky 
equation for 4 He versus cp. In all cases ce was chosen 
according to Fig. |3J The binding energy of 4 He agrees 
with the experimental one of —28.3 MeV for r\ = 1. In 
the considered range for cd , we find a unique solution for 
the cutoff choices up to A = 1.9 fm - . For A = 2.5 fm~ , 
the relation of 77 and cd is strongly non- linear and we find 
two solutions. We observed a very similar behavior, when 
the N3LO chiral interaction of IQl was augmented by 
the same 3N force. For A = 3.0 fm - , we cannot describe 
the 3 H and 4 He binding energies simultaneously. For this 
cutoff, we choose cd — 7.5, for which 77 is minimal and the 
binding energy is best described. The resulting cd/ce 
pairs are compiled in Tabled where the (*) indicates that 
the 4 He binding energy is reproduced only approximately 
as -28.8 MeV for A = 3.0 fm" 1 . 
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A [fm" 1 ] 


T 


Uow k 


:! H 
c-terms 


D-term 15-term 


T 


Mow fc 


4 He 

c-terms D-term _B-term 


1.0 


21.06 


-28.62 


0.02 


0.11 


-1.06 


38.11 


-62.18 


0.10 


0.54 


-4.87 


1.3 


25.71 


-34.14 


0.01 


1.39 


-1.46 


50.14 


-78.86 


0.19 


8.08 


-7.83 


1.6 


28.45 


-37.04 


-0.11 


0.55 


-0.32 


57.01 


-86.82 


-0.14 


3.61 


-1.94 


1.9 


30.25 


-38.66 


-0.48 


-0.50 


0.90 


60.84 


-89.50 


-1.83 


-3.48 


5.68 


2.5(a) 


33.30 


-40.94 


-2.22 


-0.11 


1.49 


67.56 


-90.97 


-11.06 


-0.41 


6.62 


2.5(6) 


33.51 


-41.29 


-2.26 


-1.42 


2.97 


68.03 


-92.86 


-11.22 


-8.67 


16.45 


3.0(*) 


36.98 


-43.91 


-4.49 


-0.73 


3.67 


78.77 


-99.03 


-22.82 


-2.63 


16.95 



TABLE II: Expectation values of the kinetic energy (T), 2N interaction (Mowfe), 27r exchange part of the 3N force (c-terms) 
and D- and _E-term for 3 H and 4 He. All energies are in MeV. 



A very important task is to estimate the size of 3N 
forces in a systematic way. We decided to calculate the 
expectation values of the 2N and the different parts of 
the 3N interactions and compare their magnitude. The 
results are summarized in Table [n] As a worst case sce- 
nario, we compare the maximum of the individual 3N 
force terms to the 2N interaction for 4 He. As expected 
from Fig. |3 for A < 2 fori. , all 3N parts are perturba- 
tive. For these cutoffs, we obtain contributions of 4—10%, 
which is comparable to 3N forces for phenomenological 
models |l3j,|22j. For larger cutoffs, the 2ir exchange con- 
tribution (c-terms) grows rapidly, which is canceled by 
the £-term. We take this as an indication that, in this 
range, our ansatz for the 3N force is not reliable. 

In summary, we have thoroughly assessed the size of 
3N forces in the V\ ow k approach. Based on the V\ ow k re- 
sults for the 3 H and 4 He binding energies, we found that 
the dependence on the cutoff is not unnaturally large 
for A ^ 1.0 fm -1 . This suggests that higher-body inter- 
actions are small. We emphasize that the large cutoff 
range, for which Vi ow fc is available, will enable similar 
studies for other low-energy observables, e.g., all binding 
and excitation energies, and that this is a powerful tool 
to isolate missing parts in effective interactions. Further- 
more, we have extended Viowfc by the leading chiral 3N 
force and fitted the two unknown parameters to the 3 H 
and 4 He binding energies. We assessed the strength of 
the 3N force by calculating expectation values of its in- 
dividual parts. By requiring that not only the sum, but 
also the individual parts are of natural size, we found 
that our ansatz for the 3N force is reliable for cutoffs 
A < 2 fm . It turned out that the 3N force contribution 
can be treated perturbatively for this range of cutoffs. 
This completes a soft nuclear interaction model, which 
will be important for many-body calculations. Applica- 
tions to symmetric nuclear matter are in preparation. 
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